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Aerosol Size Distribution Inverted from Spectral Solar Attenuation 
plus Scattering Measurements 

Inversion method_ 

The inversion technique is a least-square-method to fit measured 
optical data by optical data calculated from columnar aerosol size 
distributions via Mie theory. Starting with a single randomly se- 
lected column^ its particle concentration is chosen to fit the mea- 
sured data best. At the next step a second randomly selected co- 
lumn is added to the first column and the particle concentrations 
of both columns are chosen to fit the measured data best. In each 
of the next steps one more randomly selected column is added or 
its particle concentration is changed together with a joint change 
of all other columns to improve the best fit. The program is not 
time consuming because at each iteration step the least- square-equa- 
tion is solved analytically. The number of columns which are al- 
lowed to contribute to the size distribution^ and the size range of 
each column has to be chosen before the iteration is started. 

Test_with simulated spectral_data 

Unfortunately not only one columnar size distribution but a large 
variety of such distributions will fit the measured data. For tes- 
ting the 'uniqueness' of the inversion of spectral data to size dis- 
tributions, simulated spectral data have been used, calculated via 
Mie theory from given aerosol size distributions (including fixed re- 
fractive index) . From these tests it has been found, that the in- 
formation content of the extinction coefficients (optical depth) 
between 0.4 /am and 1.2 /im only, is insufficient with respect to 
realistic atmospheric aerosol size distributions. But if the in- 
formation content of the scattering function especially at small 
scattering angles is added, the inverted size distribution agrees 
very well with the original size distribution. Figures 1 and 2 
give for example the linear number distributions versus the par- 
ticle radius, with the straight lines being the given test size 
distributions of the simulated spectral data. They consist of ex- 
tinction coefficients at A= 0.475, 0.605, 0.743, 0.862, 1.020, 
1 . 650^*0111 and of the scattering function at the scattering angles 
2°, 6°, 10°, 14°, 18 for A= 0.605 /urn. The columns give 10 distri- 
butions resulting from 10 inversions. The 10 results are obtained 
from different random sequences of the columns during the iteration 
procedure. The figures show that the variability of the inverted 
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distributions is not 
very high, but one 
single inversion will 
not sufficiently meet 
the given size distri- 
bution. Only the ave- 
rages of the resulting 
distributions repre- 
sent both given size 
distributions (with 
and without 2 jam- 
peak) very well. The 
average of the resul- 
ting size distribu- 
tions is independent 
from. the number of in- 
versions when at least 
20 inversions are ave- 
raged. Further, the 
average is the same 
when the least-square- 
iteration is started 
with an arbitrary co- 
lumn or with any other 
distribution. This be- 
haviour of the inver- 
sion procedure justi- 
fies the application 
of the method to mea- 
sured data of spectral 
extinction coefficients 
and scattering func- 
tions . 


Sensitivity to_errors of_the_spectral data_ 

The sensitivity of the resulting aerosol size distributions to mea- 
suring errors has been tested by introducing errors into the simu- 
lated spectral data of the two test distributions presented in 
figs. 1 and 2. If the simulated data of the test distributions are 
affected with statistical errors of + 10%, the resulting size dis- 
tributions are not remarkably different from that inverted from the 
accurate data. If the errors are higher than + 20% the resulting 
size distributions are wrong with some of the columns becoming to 
high or to low what may be interpreted as holes or peaks. 

The inversion method requires knowledge of the actual value of the 
real and the imaginary part of the refractive index of the aerosol 
particles. An extensive test about the limits inbetween which an 
inaccuracy of the refractive index can be accepted is not vet made. 
But it can be said: If the true refractive index of the test distri- 
butions is 1.50 - 0.02»i and the inversion is made under the as- 
sumption, that the refractive index is 1.33 - 0 »i, the resulting 
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size distribution is wrong with fictitious holes and peaks. If the 
inversion is made under the assumption of a refractive index 
1.50 - 0»i, i.e., with only the imaginary part being incorrect, the 
resulting size distribution is more similar to the given distribu- 
tion but somewhat to inaccurate. 

Application tojmeasured data_of spectral_0£tical_deptii_and_scat-_ 
tering_f unction 


Our inversion procedure has been applied to measured data. From mid 
July to mid August 1971 ground based spectral measurements of the 
direct solar radiation and of the sky radiance at different points 

in the vicinity of 
the sun have been 
carried out in a 
savannah in the 
Southwest of Afri- 
ca near Tsumeb in 
the dry season. 

From this, reliable 
values of the aero- 
sol optical depth 
of the atmosphere 
have been derived 
at 0.465, 0.605 
0.8605 and 1.25 jam 
and values of the 
aerosol scattering 
function at the 



scatt^ring 0 angles 


10 


14 


£• » SJ f f A ** / 

18° for A=0.6o5 jum 
and X = 0.8605 pm. 
Remainig gasebus 
absorption at the 
wavelengths used 
has been taken into 
account. By means 
of a multiple-scat- 
tering-program the 

contribution of multiple scattered photons to the sky radiance has 
been taken into account, too. In figs. 3 and 4 the resulting aerosol 
size distributions of four and six consecutive days are presented. 
The exponent x at the ordinate scale has the value 5 at 28 July, 8 
at 2® July, 11 at 30 July, 14 at 31 July, 5 at 5 Aug., 8 at 6 Aug., 
11 at 8 Aug., 14 at 9 Aug., 17 at 10 Aug. and 20 at 11 Aug. The 
straight lines are power law distributions with a Junge exponent of 
3.3. They are added only to demonstrate the deviations of the real 
size distributions from the power law distribution. 


The respective types of earth surfaces which had been passed by the 
air masses in the last 3 to 4 days prior to measurement, indicated 
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on the right sides 
of the size distri- 
butions, are known 
because the trajec- 
tories of the air 
masses have been 
constructed. As to 
be seen from the 
figures, two types 
of aerosol size 
distributions have 
been found, corre- 
lated to the histo- 
ry of the air mas- 
ses. The 'arid fo- 
rest' -type has a 
peak at about 0.6 
pm particle radius 
and the 'savannah' 
type has a peak at 
about 0.9 pm par- 
ticle radius. 

Both types of aero- 
sol size distribu- 
tions are without 
any doubt bimodal 
distributions. The 
similarity of the 
different size dis- 
tributions of each 
type is - for our 
opinion - surpri- 
singly high. 
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Comparison of Extinction and Backseat tering Coefficients for Measured and 
Analytic Stratospheric Aerosol Size Distributions 

Thomas J. Swissler , Systems and Applied Sciences, Riverside, Maryland 20840 
Franklin S. Harris, Jr., Old Dominion University, Norfolk, Virginia 23508 


Theoretical Mie calculations for extinction and backseat tering coefficients 
have been performed for the analytic size distributions for stratospheric 
aerosols at 20 km described previously by Harris and Rosen (1976). The 
computations were performed under the assumption of spherical particles and 
a single refractive index characterizing the composition of the aerosols in 
the stratospheric layer. Several representative indices were chosen for 
these calculations and the wavelengths considered for examination were those 
associated with ruby and neodymium lidar systems . The ruby laser wavelength 
of 0.6942 ym and its second harmonic 0.3472 ym as well as the neodymium laser 
wavelength 1.06 ym and its second harmonic 0.53 ym were utilized. Simulations 
of aerosol extinction and backscattering at these wavelengths were performed 
for use in the interpretation of satellite measurements. The ratio of aero- 
sol extinction to backscattering may be sensitive to aerosol size distribu- 
tion or to aerosol composition. This study examines the range of sensitivi- 
ties of the theoretically calculated extinction and backscattering coeffic- 
ients for several proposed stratospheric aerosol size distributions and for 
several refractive indices. . 

The analytic models used to .describe the stratospheric aerosol size distri- 
bution have been discussed by Harris and Rosen and these models appear repre- 
sentative of experimental aerosol measurements observed during periods of 
quiescent volcanic activity. These distributions are plotted in Figure 1 in 
the form dN/dr , the aerosol number density per unit interval of radius, r. The 
particle sizes considered in this study ranged from 0.02 ym to 5.00 ym radius. 
The data on the experimental distributions presented by Harris and Rosen indi- 
cate that there are very few stratospheric aerosols outside of this size 
range. The power-law size distribution was further limited to the size range 
0.10 to 0.50 ym. Outside of these limits the predicted power-law distribution 
is not consistent with the experimental data. The total number of particles 
observed between the limits r^ and r^ is 

N ( r i’ r 2 ) = / 2 n( o dr 

Jr i 

where n(r) can also be described by 


n(r) - N q g(r) 

with the normalization constant N , consistent with the total number of particles. 
This constant has been chosen suefi that the number of particles with radii 

. r\ i . . * ~ . _ -i - '• _ 1 _ J mi „• Jj _ 


greater than 0.15 ym is equal to 1 particle cm . These size distributions are 
used with the Mie theoretical scattering terms for single particle scattering 
to define the aerosol extinction coefficient, |3. and the aerosol backscattering 
coefficient f (180°) . /* r 9 ? 

(3, - re N / r Q (2nr/A) g (r) dr 

A O^f^ e 

where Q = Mie efficiency factot for extinction 
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and 


X = wavelength (in pm) 


f (180°) = (N /k 
o 


vC 2 i(m ’ 


2nr/X 9 180°) g(r) dr 

where i(m, 2nr/A, 180°) = Mie intensity function 
k = 2 tx/A 

h* - complex refractive index. 


The Mie scattering terms in the above integrals were calculated using an 
increment in the size parameter x (where x = 2nx/X) of 0.10. Values for 
the aerosol extinction coefficients as a function of wavelength indicate 
insensitivity to variation in model size distributions with the exception 
of the aged volcanic model with its high concentration of small size 
particles. There is slight variation between models for particle refractive 
index changes. Calculated values for the backscattering coefficient as a 
function of wavelength also indicate relatively little change between particle 
size distributions with a somewhat stronger dependence on particle refractive 
index. 


If the ratio is taken of the extinction and backscattering coefficients 
a quantity is defined that is independent of the normalization used for the 
particular size distribution but still dependent on the shape of the distri- 
bution. Plots of the extinction to backscattering ratio as a function of 
wavelength are presented in Figure 2 for the aerosol composed of 75 % H^SO^ 
and 25 % E^O and in Figure 3 for aerosol having a refractive index of 1.50 
for the real part and 0.00 for the imaginary part. Better discrimination is 
seen in these ratio plots between the particle size distributions than is 
observed in calculations for either the extinction or the backscattering 
versus wavelength. For example the log normal model is distinct from the 
aged particle model, haze H and the power-law distributions. Both the haze 
H and background models have similar results because they are similar models. 


Apparently particle composition, as shown by the effect of refractive Index, 
can not easily be differentiated under any of the conditions considered in 
these calculations. However as shown in Figures 2 and 3 there is appreciable 
variation in the extinction to backscattering ratio with particle size distri- 
bution, which is also wavelength dependent. 


Reference: 

Harris, Franklin S . , Jr. and James M. Rosen, (1976): Measured and analytic 

distributions of stratospheric aerosols: a review and commentary. Atmospheric 
Aerosols , Their Optical Properties and Effects , NASA CP-2004, pp. MA2-1 to 4. 
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HYBRID METHODS ARE HELPFUL 

H.C.van de Hulst 
Leiden Netherlands 

1 . Introduction 

Methods to solve problems in radiative transfer and in multiple 
scattering exist in such a wide variety that I shall not attempt 
another review. In any practical problem the method must be chosen on 
the basis of expedience, and this in turn depends on many factors, 
such as : 

range of variables 

desired accuracy of results 

occasional or frequent computations needed 

cost, available funds 

experience and taste 

I emphasize in this talk the fact that in many situations a hybrid 
approach containing elements from different methods, though not 
"elegant", is the most practical. A normal rose or fruit tree con- 
sists of different varieties skillfully grown together because the 
desired properties of roots and fruits (or flowers) are not met in 
a single variety. The following examples may illustrate such hybrid 
methods. 

2. The X- and Y-functions 

Hybrid approaches already are helpful in the theoretical deri- 
vation, i. e. , in the process that leads to a solution in analytical 
form. The well-known X- and Y-functions for isotropic scattering 
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form a good illustration. They were introduced by Ambartsumian in 
the early forties and extensively studied by Chandrasekhar about 
1945, who defined them as solutions of certain simultaneous non- 
linear integral equations. Somewhat later I rediscovered the phy- 
sical definition as what I would now call simple examples of point- 
direction gain functions and was quite pleased to find that certain 
properties of these functions can be far more easily derived from 
these physical definitions. Since that time I do not hesitate to 
use the two approaches mixed. What happens here is that we discover 
(Fig. 1) that a derivation in mathematical physics is like a suburban 
train. Getting on at an intermediate stop, or getting off at an in- 
termediate stop, is equally authentic as taking the full ride and 
the gains in rapidity or in clarity may be considerable. 



Fig." 1 In a derivation in mathematical physics access 
from physical concepts to intermediate stops, if done 
carefully, is often advisable. Likewise interpretation 
of intermediate numerical results may be equally impor- 
tant as Interpretation of final numbers. 
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Phase functions of arbitrary form, or phase matrices with po- 
larization, require a far more elaborate set of formulae, but the 
situation is basically the same : carefully preparing the access 
at intermediate stations pays off well. 

3 , Asymptotic fitting 

A good example of a hybrid computational method is the method 
of asymptotic fitting, described elsewhere in detail. The idea is 
that by the doubling method we can in 10 steps bridge the range of 
optical thickness b = 1/32, 1/16, 1/8, 8, 16, 32, but that it 

is silly to try to carry this process to b = ca , which we can never 
reach. Instead, we wish to use the known asymptotic forms of the re- 
flection and transmission functions for sufficiently large b. The pro- 
cess by which the different approach (asymptotic theory) is grafted 
onto the numerical computation (doubling) is called asymptotic fitting. 
It uses the doubling results of 3 successive steps, say b - 8, 16 and 
32, and yields all desired functions and constants for a semi-infinite 
medium. 

4. The corner of small losses 

In the problem of reflection and transmission by a homogeneous slab 
of scattering particles, I like to imagine the results displayed in a 
diagram of a, the albedo for single scattering, against b, the optical 
thickness of the slab. At one corner^ defined by a = 1 , b = CO , the re- 
flection is complete. Small losses by escape from the bottom of the slab 
occur if b is not quite o0 . These losses are not additive and their 
combined behaviour has posed nasty problems in many papers, numerical 
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and analytical alike. Yet the situation is simple. Both these losses 
have the same dependence on angle of incidence and escape because they 
both occur in deep layers. They add in a way which is universal, except 
for scale factors, and which is shown in a particular example in Figure 2. 


loss only in olmosphere 



0.3 0.2 0.1 0 



Fig. 2 Heavy curves give the combined loss ~ non- 
ref l.ected fraction of incident flux for a finite atmos- 
phere of thickness b consisting of isotropically scatter- 
ing particles with single scattering albedo a exposed to 
vertically incident radiation. The curves near the loss- 
less corner have a universal form independent of the 
phase function. 



WC1-5 


5. Polarization in high-order scattering 

The fine visual observations of polarization of the planet 
Venus at various phase angles published in Lyot's thesis in 1929 
have outlasted 40 years before the first tentative interpretation 
could be replaced by a more definite one. Lyot conjectured that 
polarization might be completely absent in all but first-order 
scattering. This guess is not correct. I took a practicing example 
based on Rayleigh scattering in a semi-infinite atmosphere and ra- 
diation incident under 60° from the normal and reflected back to 
the source (as for a planet in opposition). The degree of polari- 
zation is 0 per cent in the first order, 9 per cent in the second 
order, and surprisingly goes up even further to 13 and 14 per cent 
in the third and fourth orders before it settles to 2,2 per cent 
in very high orders. This effect also shows up as a maximum polari- 
zation of 6.1 per cent for single scattering albedo a = 0.95. The 
simplest way to collect these results is to combine data obtained 
by very different methods, 

6. Photon path distributions 

The understanding of many problems, for instance, the formation 
of planetary absorption lines, is greatly aided by a clear knowledge 
of the photon path distributions. It is well known that such distribu- 
tions may be obtained by an inverse Laplace transform from the depen- 
dence of the reflection function on a (albedo for single scattering) . 
This knowledge may be put to use in various ways : the inverse Laplace 
transform may be applied to any form in which the reflection function 
is known, analytic, numerical, or asymptotic. A systematic exploration 
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of this possibility has made it possible to combine smoothly the 
results for very low orders, where ad-hoc calculations are fast, 
with those of high orders, where an asymptotic approach is more 
appropriate* An example is shown in Figure 3. 



Fig* 3 Path-length distribution of radiation diffusely reflected 
back in vertical direction after n scattering events from a slab 
of isotropically scattering particles of optical thickness ! ex- 
posed to vertical illumination. The exact forms, known for small 
n ? rapidly approach an asymptotic theory. 


References 

These examples have been selected from 
H.C.van de Hulst, Multiple light scattering in cloud layers, Academic 
Press, New York (in preparation), 
where further references may be found. 
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n BLUE MOON", 1975 

Raymond L . Ghuan 
Brunswick Corporation 
Costa Mesa, California 


The phenomenon of the blue moon, as was observed after the eruption of Krakatoa 
in 1883 and the extensive forest fires in Canada in 1950, had long been thought 
to be related to the anomalous extinction of light by atmospheric aerosols 
dominated by those formed by these cataclysmic events. Porch, et all, concluded 
that certain physical characteristics of atmospheric aerosols, when falling in 
prescribed ranges, could lead to anomalous scattering and the "blue moon 11 phenom- 
enon. These conditions are, (1) the geometric standard deviation of a log- 
normally distributed aerosol must be <1.5, (2) the imaginary part of the 
refractive index of the aerosol must be ^0.01, and (3) the number mean radius 
of the aerosols must lie in the range 0.4 to 0.9 jum. During a large brush fire 
in Southern California in November, 1975, with a resulting smoke cloud over 
most of the air basin for three days, automobile headlights were observed to 
be bluish (rather than reddish as might be expected in a smoke-filled atmosphere), 
especially when viewed through the rear-view mirror in a car on the freeway 
(suggesting a forward-scattering mechanism at play). With the aid of a 10-stage 
cascade impact or. capable of generating real-time in-situ data on aerosol mass 
distribution by size from 0.05 to 25 jum, a number of aerosol size distributions 
were measured from a car on the freeway, yielding some interesting results^. 

The dominant feature of these distributions, obtained repeatedly over a two-day 
period, was a tri-modal distribution with a middle mode, centered mass-wise 
about 0.9 -urn diameter, distinctly separated from two other modes with mass-median 
diameters of 0.12 Aim and 4.3 jum. Figure 1 shows one such distribution, shown as 
differential normalized mass concentration versus aerosol diameter. By treating 
each of the three as a separate mode, as indicated by the dashed portions of the 
curve, three distinct log-normal distributions can be discerned, as shown in 
Figure 2. From these the relevant aerosol characteristics are derived and shown 
in the table below. 


Mode 

dm 

(um) 

°g 

dg 

d 

dm 

da 

A 


C 

(g/m 3 ) 

1 

0.12 

1.75 

0.047 

0.055 

0,075 

0.064 

8.5 

1.08 

x 10-4 

2 

0.94 

1.40 

0.67 

0.71 

0.79 

0.75 

23.4 

5.25 

x 10-5 

3 

4.3 

1.32 

3.4 

3.5 

3.8 

3.7 

34.4 

2.41 

x 10-5 


dm: 

Mass median diameter, micrometer 

A: 

Shape parameter, (2 log^og) 

^g: 

Geometric standard deviation 

C: 

Mass concentration, g/m3 

dg: 

Geometric mean diameter 

Oi 0 : 

Size parameter , 7rdg/\ 

3: 

Count mean diameter 

X: 

Wavelength 

dm: 

Mass average diameter 

kg ; 

Extinction coefficient, m^ 

da : 

Area average diameter 

k s : 

Scattering coefficient 
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Using the computation results of Foitzik, et al , as presented in the Handbook 
on Aerosols by Dennis^*, the extinction coefficients and the forward scattering 
(10 degrees) coefficients for the three modes have been calculated, for three 
wavelengths, and the results are shown below: 


Mode 


d 

A 


A 

k e 


k s 


1 

1 . 75 

0.055 

8.5 

0.33 

450 

~0 


~o 






0.27 

550 

~ 0 


~0 






0.23 

650 

~0 


~0 


2 

1.40 

0.71 

23.4 

4.67 

450 

2 .3 x 

10-4 

3.0 x 

10“ 5 





3.82 

550 

2.4 x 

10" 2 3 4 

2.4 x 

10-5 





3.23 

650 

2 .4 x 

10“ 4 

2 . 0 x 

10-5 

3 

1.32 

3.5 

34.4 

23.8 

450 

1.4 x 

10-5 

9 .4 x 

10-6 





19.5 

550 

1.4 x 

10-5 

7.0 x 

10-6 





16.5 

650 

1.4 x 

10-5 

5 . 8 x 

10-6 


It is seen that Mode 1 is largely below the Mie scattering regime. Mode 2 
satisfies the conditions derived by Porch^. Since forward scattering can be 
seen to be a dominant mechanism, the imaginary part of the refractive index 
has little influence. There is a decrease in the extinction coefficient with 
decreasing wavelength; but the main influence in the "blue headlight" observa- 
tion appears to be the enhanced forward scattering with decreasing wavelength. 


REFERENCES 

1. William M. Porch, David S. Ensor, Robert J. Char Is on, and Jost Heintzenberg, 
Applied Optics, 12, pp, 34-36, Jan. 1973. 

2. Raymond L. Chuan, "Real Time Measurement of Aerosol Mass Distribution By Size" 
Ninth Aerosol Technology Meeting, October 21-22, 1976, Battelle Columbus 
Laboratories . 

3. I. Foitzik, G. Hebermehl, and D. Spankuch, Optik, 23, 272 (1965/66). 

4. Richard Dennis (Ed), "Handbook on Aerosols", Technical Information Center, 
Office of Public Affairs, U.S.E.R.D.A. , TID -26608 (NTIS) , 1976. 
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FIGURE 2. PARTICLE LOG-NORMAL DISTRIBUTION 



PD2-1 

Absorption Refractive Index of Aerosols in the 
Denver Pollution Cloud in the Interval .53— .69 \X 


J. J. DeLuisi 
NOAA/ ERL-AJRL , Boulder, CO 
J\ E% Bonelli and C. E % Shelden 
NCAR, Boulder , CO 


OVERVIEW 

Wideband spectral measurements of the global and direct solar flux at the 
surface were obtained at the Denver Trout Farm in conjunction with the EPA 
Denver Brown Cloud experiment on 20-21, November 1973.. These measurements were 
analyzed by a simple variational technique to yield the increase in the rate 
of absorption and scattering by the Denver pollution cloud when the observed 
optical depth increased by approximately 0.06 . The imaginary term of the com- 
plex refractive index of the aerosols is deduced from theoretical calculation 
data. 

PROCEDURE 

On 20 and 21 November 1973 , observations of the direct and global solar 
flux at the earth* s surface were obtained near the Denver Trout Farm. The di- 
rect beam is measured normal to the solar beam. The global flux is the sum 
of the skylight radiation and the direct beam incident on a horizontal re- 
ceiver. 

During the two-day period of measurement the optical depth of the pollu- 
tion layer increased by approximately 0. 06 from the first to the second day- 
and the solar flux reaching the earth 9 s surface decreased by approximately 
4.8 % over the same time interval* The measured decrease in flux was ascribed 
to increased scattering and absorption by aerosols and gases in the polluted 
boundary layer. 

The solar global flux impinging on a flat-plate receiver at the earth 9 s 
surface is given by P where 

P( y 0 , rn, T, r) = D ( V 0 , T) + S < P Q , m, T, r) (1) 


where D is the direct radiation given by the expression 

~ x/ u 

D C^O »' T > = V 0 F 0 e 0 * 


( 2 ) 


and where 

Vq - cosine of the solar zenith angle ; F ^ = the incident solar flux 
m - complex refractive index of aerosols ; r = surface reflectivity , 
T - total optical depth of the atmospheric medium . 


and 
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The S term in Eq (1) stands for the flux of the diffusely transmitted radia- 
tion which in our approximation is radiation scattered once downward after 
having been scattered from the direct beam , and back- scattered once downward 
after having been reflected by the earth's surface . Therefore , S depends on 
(1) surface reflectivity , (2) scattering and absorption optical depth of aero- 

sols and gases , (3) solar zenith angle, and (4) the solar flux F incident at 
the top of the optical medium . Surface reflectivity is important when there 
is substantial snow cover as in this work ( r^ . 5 ) 

The S term can be separated into S and S which are the forward scattered 
flux removed from the direct beam and the downward back-scattered flux (down- 
welling) arising from the flux reflected once by the surface , respectively . 

If we take the derivative of Eg (1) with respect to , the optical depth 
of aerosols we have 

iS „i»_ + iS£ + lSt. <3 > 

St a 3r a 3T a 3 - T a 

We can express Eq (2J by an approximation having the form 

P = y 0 F 0 e" X /y 0 + aK)i 0 F 0 (l-e“ T 1 % + r(l-a)‘ KP q F o (l-e~ T 1 % (4) 

where a=0.82 (at a wavelength of 0.525 V ) is the fraction of power scattered 
downward and K, the unknown variable , is the fraction of power removed from the 
incident flux that is available for scattering. The term a is quite compli- 
cated for aerosol scattering, but it can be estimated in a formal way using 
Mie theory and a known aerosol size distribution. 

Differentiating Eq (4) according to Eq (3), and solving for K we have 


K = 


1 

a+(l-a)r 



(5) 


where -AD = F Q e T ^ V °Ax d (for a horizontal receiver) 

which is the change in power removed from the direct beam in finite difference 
form . 

MEASUREMENTS 

Radiant flux measurements were made with commercially available radiometers . 

The total flux was measured using three Eppley Model 2 pyranometers , one in each 
spectral band of interest . The direct normal flux was measured with an Eppley 
model 15 pyrhel iometer equipped with a three position rotary filter holder . WG7 ,OGl 
and RG8 Schott filters provided the spectral resolution . The pyranometers and 
pyrheliometer were not intercalibrated , nor do they need to be . 
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Figure 1 provides a summary of the analysis , following the procedure 
described above . It is evident from the figure that comparatively strong 
absorption is taking place in the OG1 and RG8 bands . 

RESULTS 

The work by Herman et al, (1975) provided us with a means to estimate a 
refractive index for the particulate materials in the brown cloud , We assumed 
all absorption in the OG1-RG8 band is due to particles having a size distri- 
bution approximated by dN/dlogr ~ bl where N is particle number, 1 is par- 
ticle diameter , and b is a constant . Assuming also that the real term of the 
complex refractive index is 1,5 (Bhardwaja et al 9 , 1974), then the imaginary 
term can be deduced from functions similar to those of Herman et al . , (1975): 
percent absorption incident solar power versus absorption refractive index for 
a specified solar zenith angle, surface albedo and aerosol optical depth . Fol- 
lowing this procedure , we arrived at an approximate value of -0.03 for the imag- 
inary term. This value is rather large, but is consistent with the strong 
absorption occurring within the cloud . Simultaneous radiation absorption by 
molecules will fictitiously increase the magnitude of the absorption refractive 
index. We believe it is impractical to estimate a figure of uncertainty on the 
value (in view of the uncertainties connected with the data) , other 
than to state that it may be quite large, and that our estimate is probably 
an upper limit. 

Our refractive index can be compared with the results of other investiga- 
tions of an urban atmosphere* Eiden (1966) deduced a range of -0.01, -0.1 for 
Mainz from polarization data, and Lin et al . , (1973) deduced values ranging 
from - 0,028 to -0..05 for New York City air, using the opal glass technique. 
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A TWO-DIMENSIONAL STRATOSPHERIC MODEL OF THE DISPERSION OF AEROSOLS 

FROM THE FUEGO VOLCANIC ERUPTION 

Ellis E. Remsberg, NASA Langley Research Center, Hampton, Virginia 23665 
Carolyn F. Jones, Old Dominion University, Norfolk, Virginia 23508 
Jae H. Park, College of William and Mary, Williamsburg , Virginia 23185 


Observational data of the pre- and post-volcanic aerosols have been used 
in conjunction with predictions of a 2-D circulation model to gain better 
understanding of the transport, chemical, and sedimentation processes for 
the stratospheric layer. This paper specifically concerns the Fuego 
eruptions in Guatemala (15°N) in October 197*+ • A lidar measure of 
aerosol mixing ratio used for observational data in this paper is the 
scattering ratio, R s = 1 + f a / f m , where f a and % are the aerosol 
and molecular backscattering functions, respectively. The f values 
are products of the species cross section and number density. Figure 1 
shows two R s profiles obtained by lidar at Hampton, Virginia. The plot 
for January 2, 1975, represents an enhanced aerosol layer due to Fuego 
volcanic activity. The February 19, 1976, profile resembles a near- 
background aerosol level and shows the depletion of the January 2 layer 
over 13 months (ref. l). Any R s value greater than one represents 
backscattering from aerosols. If the aerosol cross section is constant 
with height , the scattering ratio profile is a direct measure of aerosol 
number density. 

The aerosol distribution is predicted by the model in the form of 
y(t) = + (dx/dt)At where the rate of change of aerosol mixing 

ratio is defined by dx/dt = Ox/3t) c h + Ox/3t ) gr + Ox/3t) tr + Ox/9t) sed 
X is the aerosol mixing ratio, (3x/3' t ) ch is the chemical term, 

(3x/3t)g r is the aerosol growth term, Ox/3t)-fc r is the transport term, 
and (9x/9t) sed is the sedimentation term. The numerical integration 
is performed from initial time, t 0 , to final time, t, at intervals 
of At. Since the gas_to_aerosol chemistry is not well understood, the 
modeling attempt in this paper does not consider the term (9x/9t ) c h* 

Also, because the processes of aerosol growth are only understood 
qualitatively (ref. 2), they have not been incorporated into the model 
calculations at this time. 

The dynamical model developed by Louis (ref. 3) includes adve.ct.ion by 
the mean meridional circulation and diffusion by large-scale eddies. 

The model extends from 0 to 50 km in altitude with a grid spacing of 
1 km, and from 90°N to 90°S with a grid spacing of 5° • For simulations 
presented in this paper, the circulation has been specified by monthly 
mean winds and eddy diffusion parameters derived from Louis’ Model II. 

Mean winds for the winter season are shown in figure 2. Louis’ model 
has successfully approximated the distribution of radioactive debris in 
the stratosphere. In particular, the analysis of a volcanic event 


*Jae Park is supported by NSF DES75-l*+228 for this work. 
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represents dispersion from a point source, similar to that for radio- 
active bomb debris. 

Hunten (ref. 4 ) has discussed the importance of aerosol sedimentation 
rates for determining the residence times of volcanic aerosol layers. 

Aerosol fall speeds for various aerosol sizes tabulated by Kasten (ref. 5) 
for particle densities of 1.5 g/crn-^ have been applied to an initial aero- 
sol size distribution. An estimate of the Fuego aerosol size distribution 
has been adopted from the 1963 measurements by Mossop (ref. 6) which were 
obtained 1 month after the eruption of Mount Agung. The total size dis- 
tribution has been divided into four size ranges with mean radii of 
0.l6 pm, 0.32 pm, 0.53 pm, and 0,93 pm. Distributions of mass mixing 
ratios of four different sizes are then calculated as functions of alti- 
tude and time. 

Initially, it is assumed that any substantial deviations between the model 

results and the observed aerosol profiles can be accounted for within 

the uncertainties in the aerosol sedimentation, growth rates, and chemistry. 

High resolution infrared satellite photographs (released by the national 
Oceanic and Atmospheric Administration in Rockville, Maryland) were 
employed to estimate the initial location (shown as a dot in fig. 2), 
size, and direction of the dust clouds for several days after the eruptions 
of October lit and 17, 197^* The model is started on October 19 at 1200Z 
to allow for some spread of the cloud.* Data taken hy lidar .at Hawaii 
(ref. 7) on October 29 were applied to verify the initial vertical profile 
of the dust layer. The vertical width of the layer at half maximum as 
computed by the model for October 29 at 20°W was compared with the observa- 
tions (half-width of 0.8 km) at Mauna Loa Observation in Hawaii. The 
shapes of the model profiles at all latitudes where aerosol had been 
transported in that 10-day period were adjusted to agree with the Hawaii 
observations. These adjusted profiles then represent the initial con- 
ditions for the aerosol source. Amounts of injected material were deter- 
mined from estimates by'Cadle, et al. (ref. 8). A constant mass mixing 
ratio of 2 x 10-10 was assumed for the model at a lower boundary of 10 km. 
Dustsonde data from Wyoming (ref. 9) support the assumption that the 
aerosol mixing ratio in the upper troposphere remains almost constant 
with time after the eruption. 

Figure 3 displays the integrated aerosol mass density between l6 and 21 km 
as a function of time after the eruption. The solid line represents the 
lidar data from Hampton, Virginia (37°N) • The lidar data were obtained 
by summations over 1 km increments of (R S (Z) - l) W m (Z) where R S (Z) 
is the average scattering ratio and N m (Z) is the molecular number 
density. The quantity on the ordinate is then equivalent to (o & /o m ) N a (Z) 
where a a and a m are aerosol and molecular cross sections, respectively, 
and N a (Z) is the aerosol number density which is directly proportional 

to aerosol mass density. Wo attempt has been made to actually compute 
mass densities from the lidar data. The model results in figure 3 are 
an average for 35°W and U0°N latitude. The model quantity is a little 
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different from observed quantities because it does not contain the addi- 
tional effects of a time-varying size distribution as does the scattering 
ratio lidar measurements which are more sensitive to larger size particles. 
The circles in the figure are for a log-normal background type aerosol 
size distribution (ref. 10) where the mean radius is 0.0725 pm. The 
triangles are for a volcanic aerosol size distribution taken from Mo s sop 
(ref. 6). The model results in figure 3 clearly show that vertical aero- 
sol profiles and aerosol residence times are sensitive to the sedimentation 
term, as well as the transport term. 

The peak magnitudes of the lidar and model data have been arbitrarily 
adjusted to afford a better comparison of the time rate of change of the 
aerosol column load. The time of the occurrence of the maximum aerosol 
load has been simulated very well by the model, indicating that the 
meridional transport for the first few months is reasonable. The l/e 
decay time for the integrated aerosol column density from lidar measure- 
ments (16-21 km) after February 1975 is about 10 months. It is important 
to note that longitudinal movements of aerosol clouds, displayed in the 
observational data as periodic peaks, are not simulated in the present 
zonally averaged model. 

Figure h shows the initial size distribution adopted from Mossop for the 
Agung eruption for 20 km and 15°N to 35°S latitude and then compares sub- 
sequent measured size distributions with those determined from the model by 
including sedimentation processes. In general, the comparison is reason- 
able. The initial size distribution employed in this model may, however, 
contain too many large particles. The sedimentation rates may be too 
rapid, which, would affect the results shown in figure 3 as well. Compara- 
tive vertical profiles of mass mixing ratios will not be presented in 
this paper. 

The initial study into the dispersion of volcanic particles has shown that 
transport alone does not account for the spread to other latitudes and 
altitudes. To obtain more accurate sedimentation rates, it will be 
necessary to have more post-volcanic- particle size distribution data. 

The aerosol chemistry may be of importance and should be included in 
further studies. Quantitative knowledge of aerosol growth by gas con- 
densation must also be considered. 
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ON THE VARIATION OF THE AEROSOL VOLUME 
TO LIGHT SCATTERING COEFFICIENT 

by 

George M. Sverdrup* and Kenneth T. Whitby 

Department of Meehan i cal Engineer!’ ng 
University of Minnesota 

*Present Address: Battel le, Columbus Laboratories 

Columbus, Ohio 


A number of investigations have been conducted relating aerosol 
mass concentration to the aerosol light scattering coefficient b S p (m"'). 
These studies, reviewed by Chari son (1969), suggested, that the ratio of 
aerosol mass or volume to bsp was approximately a constant, and that this 
relationship could be used on a practical basis to infer the aerosol mass 
concentration from b S p. 

The present work examines the variation of the ratio, Kp, of 
submicron aerosol volume, V3 - (ym^/cmS) ,. to b s p. The experimental data 
were acquired in the California Aerosol Characterization Experiment (ACHEX) 
during 1972 and 1973. V3 - was determined using the Minnesota Aerosol 
Analyzing System (Whitby, et al . , 1975). An integrating nephelometer was 
used to determine b S p. 

The variation of K p as a function of the total aerosol number 
concentration, CNC (No./cm3), determined with a condensation nuclei counter 
is shown in Figures 1 and 2. K p ranges from 5 to 80. Near combustion 
sources large quantities of particles with sizes below the light scattering 
range were added to the distribution increasing V3 - significantly but not 
affecting ;'i>sp appreciably. Thus as CNC rose, Kp became greater. 

the data fall into groups which can be identified with general 
location. Six categories are shown in Figures 1 and 2. Aged aerosols have 
values of Kp in the range 6 to 20 ym^/cm^ x m x 10"4. Fresher aerosols, 
e.g., categories D and F, have values of Kp in the range 20 to 80. The 
two nonurban categories have values of CNC less than 20,000 cm"3. Urban 
areas not influenced by strong combustion sources have concentrations 
between 20,000 and 250,000 cm" 3. 

The variation of K p was examined theoretically by modeling the 
aerosol size distribution with a composite distribution made up of 3 log- 
normal functions. These components correspond to the nuclei, accumulation 
and coarse particle size distribution modes (Whitby, 1975). The variation 
of Kp was studied by computing the value of Kp from the model distribution 
using Mie theory to predict b S p. A two-dimensional nephelometer geometry 
was used. The light scattering efficiency was computed by integration over 
the nominal spectral response of the nephelometer using a refractive index 
of 1.50-0.02 i. The initial model distribution was the fit to the grand 
average measured distribution from Goldstone, CA. The concentration of the 
nuclei mode was then varied in order to simulate different aerosol sources. 



PDk-2 


The resulting variation of Kp is shown as the dashed curve in Figure 1. 

The conclusions to be drawn from this work are 

(1) The large variations in Kp (i.e., 5 to 80) preclude 
its use as a good indicator of submicron aerosol 
volume or mass concentration. 

(2) Using the particle number concentration and Kp, 
aerosols can be classified according to general 
location. 

(3) Variations in aerosol size distributions are enough 
to account for the observed variation of Kp. 
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Spectrophone In Situ Measurements of the 
Absorption of Visible Light by Aerosols 

R. W. Terhune and J. E. Anderson 
Research Staff, Ford Motor Company, Dearborn, Michigan 48121 

In extinction measurements of visible light by atmospheric aerosols, the 
effects of absorption are usually masked by scattering. Separate absorption 
measurements would aid greatly in interpreting extinction data and character- 
izing aerosols. It would also allow one to evaluate the contribution of vis- 
ible sunlight absorption by atmospheric aerosols to the radiation balance of 
12 

the earth. * Aerosol absorption properties are usually deduced from measure- 
ments of the optical properties of collected samples and measured particle 

3 

size distributions. In a paper at this meeting, Patterson, Gillette and 
Stockton report on their recent measurements of the complex index of refrac- 
tion of aerosols from 300 to 700 nm. Also, Bruce and Pinnick^ report at this 
meeting some infrared spectrophone in - situ aerosol absorption measurements. 

We report here in situ measurements of the absorption by atmospheric 
aerosols in the visible near 514.5 nm. We used a 1/2 watt 514.5 nm argon 
laser beam chopped at 4000 Hz to drive a spectrophone. Light absorption led 
to localized heating. The associated gas expansion in the cavity generated 
a sound wave which was detected with a microphone. 

In the case of absorption of light by aerosols, sound generation is 
limited by the rate at which heat is transferred to the gas. The thermal 

time constant of the particles is proportional to the square of their diam- 
6 

eter. At 4000 Hz, heat transfer is efficient for spherical particles less 
than about 3 microns in diameter. For atmospheric measurements, it would be 
desirable to detect the absorption in particles up to about 20 microns by 
operating at 90 Hz or less. One can use this particle-size dependence of the 
time constant to determine the absorption contributions from different par- 


ticles. 



PD5-2 


In order to obtain enhanced sensitivity, a resonant acoustic cavity was 

7 8 9 

used in the spectrophone. ’ * A cylindrical brass cavity* 5-cm in diameter 
and 15-cm long* was used in its lowest order circumferential mode which had a 
Q of 250. The laser beam was passed through the cavity near one wall, enter- 
ing and leaving through 8-mm diameter holes in the end walls. A 1-inch B&K 
condenser microphone with a sensitivity of 42 mV/Pa was mounted flush with the 
cavity wall on the opposite side. Synchronous detection was used. An acous- 
tic signal of about 1.1 Pascals rms per m ^ absorption per watt of laser power 

- 1/2 

was obtained. An rms noise level of about 0.01 pV sec was observed after 
synchronous detection. Acoustic isolation was achieved by placing the cav- 
ity inside three sound insulating boxes, one inside the other. Measurements 
were made with gas flowing at a rate of 2 Jl/min. The system was calibrated 

by observing the absorption of 30 ppm NO^ in air which was measured to be 
-4 -1 

around 180 x 10 m . The above numbers indicate a. system detection limit 

-7 -1 

(S/N ratio = 1) of 2 x 10 m in one second with a single pass of a one-watt 
laser. ^ The acoustic circuit including the detection system design has not 
been optimized. 

Measurements with an integrating nephelometer indicated an extinction of 

around 2.5 times Rayleigh scattering in the room air . The absorption measure- 

-5 —1 

ments in the room air shown in the figure indicate an absorption of 10 m , 
probably due to aerosols with diameters less than 3 microns. The residual 
absorption after filtering is about what one would expect from 3 ppb of 

A quick series of experiments were carried out with NaCl and NiCl^ aero- 
sols made by aspirating water solutions. Other experiments used cigarette 
smoke. The aerosols were blown through the cavity. Scattering measurements 
were made with an integrating nephelometer sampling the exiting stream of 
aerosols. The results are summarized in the table. 
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Recorder traces showing the observed synchronously detected spectrophone 
signal with and without aerosols, A millipore filter with a 10-nm pore size 
was used. The large fluctuations in the signals with ulfi.lt.ered air can 
probably be associated with fluctuations in the number of particles illumi- 
nated by the 1-mm diameter laser beam in the spectrophone. 0.12 y volts 

corresponds to an extinction coefficient due to absorption of 0.5 x 10 ^ 

-2 

m . The extinction coefficient due to Raleigh scattering in air at this 
wavelength is about 10 m **". 
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Aerosol Extinction Coefficients Relative to Rayleigh Scattering 




Scattering 

Absorption 

Ratio 

NaCl 


50 

<0.2 

<0.004 

NiCl 2 


50 

1.5 

0.03 

Cigarette Smoke 

through 
off tip 

250 

15 

200 

100 

0.8 

7 


The NaCl experiments demonstrate the insensitivity of the spectrophone to 

scattering per se . Here only transparent aerosols are expected. Our analysis 

— 

indicates that only about 10 of the scattered light would be effective in 
generating acoustic signals following absorption in the cavity walls. In 
the case of NiCJ^, a small amount of absorption was found as expected* 

For the aerosol drawn through the cigarette, the absorption and scatter- 
ing were observed to be about equal, which one would expect for an aerosol of 
medium to large particles of optically lossy material. For the aerosol off 
the tip of the cigarette, the absorption is much larger than the scattering 
which is characteristic of small particles of optically lossy materials. 

These results are generally consistent with the differences in the particle 
size distributions observed by Whitby el al. f or cigarette aerosols in the 

two cases . 
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THE KRAKATOA VOLCANIC TURBIDITY - 

A REASSESSMENT 

Frederic E. Volz 

Atmospheric Optics Branch, Air Force Geophysics Laboratory 
L. G. Hanscom AFB, Bedford, MA 01731 

The explosions of Krakatoa volcano in August 1883 (Symons 1888) are considered 
to have been the strongest volcanic event in the last 150 years. Tremendous 
amounts of ash reduced the solar radiation for some years. The general effect 
of volcanic dust on climate is still unsettled. Nevertheless, using dust 
data since the Krakatoa eruption, Oliver (1976) could obtain close agreement 
with observed northern hemispheric surface temperature variations. 

A great number of worldwide reports of visual observations of abnormal twi- 
light and skylight phenomena due to the Krakatoa dust has been summarized 
and discussed (Symons 1888, and others). There have been recent attempts to 
quantitatively discuss some aspect of the dust event (De i rmendj ian 1973, Dyer 
1 97*0 but measurements hitherto not considered make possible a more detailed 
and balanced account. Also, more recent data allow to develop more generally 
the fate of stratospheric dust of a tropical source. 

In the period of the Krakatoa event (since 1875) solar radiation was measured 
only at Montpellier (42N), at the French Mediterranean. Monthly averages of 
the noontime intensities presented by Kimball (1918) and used again by Dyer 
(1 97*0 show the greatest turbidity about 2 years after the explosion. Having 
been intrigued by the long delay and poss i bi 1 i ty of large seasonal variations, 

I examined the original data and came to the following conclusions: (1) Turbidity 

at this sea level locality varied greatly so that averages of months with few 
measurements are very questionable. (2) As shown in Fig. la, volcanic aerosol 
optical thickness was about 0.15 during 1884 and 0.10 in 1886 when other 
sky phenomena were said to have faded out (see Fig. 1c and d), and seem to 
have been recognizable as late as 1888 and 1889 (eruptions by other volcanoes 
in 1886 and 1888 seem to have been minor). Comparing radiation maxima only, 
volcanic turbidity appears to have been about 30% larger (Fig. lb). 

Some measurements of the position of neutral points of sky polarization in 
Germany (Busch 1 889 ) between late 1886 and late 1888 (Figs. Id and 2) confirm 
that volcanic turbidity ( T ^ « 0 . 02 ) still was present by early 1888, Also, 
the moon was still very dark in eclipses of 1888 and 1 889 - 

Repo r ts from the inner tropics of heavy sky haze, dim and colored suns even 
at noontime, and the Bishop ring have, as well known, shown that the fresh 
Krakatoa dust cloud circled the earth twice from east to west in 2 weeks 
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time. The few reports of the sun disappearing before setting in striated 
haze of strange colors lead, based on our corresponding observations in 
heavy summer haze, to estimates of 0.45, 0. 60 and 1.6 (De i rmendj ian 
1973 used 0.65). However, these seem to have been extreme cases before the 
dust got more evenly distributed. Indeed, 2 months after the eruption there 
weren't any more relevant reports from the tropics, although Bishop rings 
and purple lights lasted at least until 1886. 

We, therefore, evaluated as in our Katmai study for higher latitude sources 
(Volz, 1975) nuclear debris data from sources in the near-tropical strato- 
sphere. -As Fig, 3 shows the steady state vertical debris load increases 
toward ‘the poles by more than a factor 2. This apparent paradox follows 
from the fact that the peak mixing ratio in the relatively thin debris layer, 
which slants from about 2 k km in the tropics to \k km in the Arctic, decreases 
along this profile less than the air density. Indeed, the evaluation of solar 
data of the Agung period by Dyer and Hicks (1 965) Suggests constant average 
turbidity between 18N and 81 N, and shows notably higher .turbidity at the 
South Pole than at the same time at 30-kOS (the high turbidity derived for 
0-10S is, at least with regard to the Congo data, very questionable as are 
the northern hemisphere minima). 

The decrease of midlatitude (~ k^N) turbidity from great tropical volcanic 
eruptions is shown in Fig. k, superimposed by the year of eruption, and taking 
data by Dyer (197*0 into account. About 1.5 years after the eruption, 

Krakatoa turbidity is higher by a factor 5 than Agung turbidity (crude global 
average), and then abates with a residence time of about 1.5 years. That, 
according to Tidar data, the small Fuego turbidity vanished much faster is 
attributable to its 1 confinement to altitudes below 19-20 km even. in the 
tropics. Tropical twilight observations gave altitudes of 25-35 km for the 
Krakatoa cloud, and indeed easterly winds of 30m/sec as derived from the drift 
of the fresh cloud do not seem to normally exist below about 2k km (Belmont 
and Dartt, 1 970) . 

There are more problem areas of the optics of the Krakatoa dust, which briefly 
will be mentioned. Thus, the influence of heavy sky haze on twilight phenomena 
notably the possibility to derive the layer altitude from visual observations 
when extinction becomes dominant, deserves attention. Very interesting is also 
the size distribution. While the radius derived from the globally reported 
Bishop ring is about 0.6pm, this size cannot - as Deirmendjian (1973) pointed 
out - have been very predominant because the sun was normal red, and not green, 
at setting except early in the tropics, Otherwise, this Bishop ring might be 
related to the di fference of the Arago neutral points in Krakatoa and Katmai 
twi 1 i ghts ( F i g. 2) . 

















pig. t . Comnarison of turbidity at about 45N from tropical 
stratospheric volcanic dust source. 

Global aerosol load* 0»3 .years after eruption, for 
Krakatoa 34, Katmai 13, Agu ng 9, Fuego 1.5(mio tons). 
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Multiple Scattering in the Sky Radiance 


by Bary de, E. and Bui 1 rich , K. 

Institut fur Meteorol ogie, Johannes Gutenberg-Uni vers i tat, 
65oo Mainz, Germany 


Multiple scattering is defined as scattering of higher 
order plus reflected radiation from the surface within 
the solar wave length range. 

Distinct ways of the solution of the equation of radiative 
transfer, containing the multiple scattering, make it possible 
to calculate the sky radiance and its degree of polarization 
at any points in the sky for any altitudes in the atmosphere. 
That means under distinct assumptions based on measurements 
of the optical properties of -the particulate matter in the 
atmosphere, it is possible to compute the upward and down- 
ward radiative fluxes and their divergences. These diver- 
gences can be expressed as heating or cooling rates respectively 
at different heights in the atmosphere. These results are 
of great interest for the energy budget of the atmosphere. 

The following pictures give examples a) of the amounts of the 
multiple scattering or the scattering of higher order re- 
spectively partly in the form relative to the amount of the 
primary scattering, b) They show the distribution of the 
portion of the multiple scattering in the sky. c) They de- 
monstrate the influence of distinct optical properties of 
the particulate matter, namely the aerosol particles. 
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Fig. 1 Portions of the first, second and third scattering 
processes of sky light. 

Fig. 2 Ratio of sky light including scattering of higher 
order over primary scattering only. T turbidity 
factor, m complex refractive index, wavelength 
A= 0.55 microns, surface albedo A = o, zenith- 
distance of the sun = 37°. The computated results 
are valid for sun_s vertical. Eschelbach, G. 1974. 

Fig. 3 Percentage ratio of primary PS to total scattering TS, 
as a function of the scattering angle and the 
wavelength /V at the solar elevation of 24°. Surface 
albedo A = o.o. No aerosol absorption, the hatched 
surface marks the difference in the values of the sun's 
vertical and the solar almucantar. Turbidity factor 
T = 2. Bary de, E. and Eschelbach, G. 1974. 


Fig. ’Three dimensional plot of surface of equal percentage of 
the ratio PS/TS. The orthogonal coordinate system has 
the angular distance from the sun ( ) as x-axis, the 

wavelength (A) as y-axis and the turbidity factor (T) 
as z-axis. S.V. means sun's vertical . Surface albedo 
A = o.25, complex refractive index m = 1 . 5-o.oi 
Bary de, E. 1973. 
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